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Abstract. We study the one-point probability distribution functions (PDFs) of the 
pecuHar velocity and the density fluctuation in a cosmological fluid. Within the 
perturbative approach to the structure formation scenario, the effect of "pressure" has 
recently been an area of research interest. The velocity dispersion of the cosmological 
fluid creates effective "pressure" or viscosity terms. From this viewpoint, because the 
pressure reflects a nonlinear effect of the motion of the fluid, the pressure model would 
include nonlinear effects. Here we analyze the Lagrangian linear perturbation PDFs 
for both the Zel'dovich approximation and the pressure model. We find that the PDFs 
of the peculiar velocity remain Gaussian, even if we consider the pressure. For the 
PDFs of the density fluctuation, the occurrence of non-Gaussianity depends on the 
"equation of state" for the fluid. Therefore we distinguish the "equation of state" 
using the PDFs. 
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1. Introduction 

The scenario for tlie formation of large-scale structure in the Universe is based on the 
gravitational instability of the primordial density fluctuation. This fluctuation may have 
been generated from the quantum fluctuations during the inflation phase [H El E! • The 
fluctuation is spontaneously grown by its self-gravitational instability. To describe the 
evolution of the density fluctuation, many approaches have been carried out. 

For a perturbative approach, the Lagrangian description provides a relatively 
accurate model even in a quasi-linear stage. Zel'dovich |3] originally proposed a linear 
Lagrangian approximation for dust (pressureless) fluid. This approximation is called 
the Zel'dovich approximation (ZA) [IliaillElinilZllHll3linilIIl- ZA describes the 
evolution of density fluctuation better than the Eulerian approximation ^21 ^1 El- 
Although ZA gives an accurate description until a quasi-linear regime develops, ZA 
cannot describe the model after the formation of caustics. In ZA, even after the 
formation of caustics, the fluid elements keep moving in the direction set up by the initial 
condition. In order to proceed with a hydrodynamical description without the formation 
of caustics, although modified models, such as the "adhesion approximation" jTH] and 
the "Truncated Zel'dovich approximation" [ini Ej were proposed, the physical origin 
of the modification has not yet been clarified, the physical origin of the modification has 
not yet been clarified. 

Recently, instead of the dust fluid, pressure-supported fluid has been considered. 
The coUisionless Boltzmann equation PH] describes the motion of matter in phase space. 
The basic equations of hydrodynamics are obtained by integrating the coUisionless 
Boltzmann equation over velocity space. In past approximations, such as ZA and its 
modified models, velocity dispersion was ignored. Buchert and Dominguez [19J argued 
that the effect of the velocity dispersion should be noticeable beyond the caustics. They 
showed that when the velocity dispersion can be considered isotropic, it gives effective 
"pressure" or viscosity terms. Furthermore, they argued for a relation between mass 
density p and pressure P, i.e., an "equation of state." If the relation between the density 
of matter and pressure can be regarded as barotropic, they showed that the "equation of 
state" should take the form P oc p^/^. Buchert et al. [201 showed how the viscosity term 
is generated by the effective pressure of a fluid under the assumption that the peculiar 
acceleration is parallel to the peculiar velocity. Dominguez |211 introduced the idea 
that the evolution equations for the matter fields are smoothed over a smoothing length; 
then the viscosity term in the "adhesion approximation" can be derived by the expansion 
of coarse-grained equations. Recently, Buchert and Dominguez j2S] discussed the origin 
of the viscosity term and the extension of the Lagrangian perturbation theory. 

Departing from these points of view, the extension of the Lagrangian perturbation 
theory to cosmological fluids with pressure has been considered. Adler and Buchert [21] 
formulated the Lagrangian perturbation theory for a barotropic fluid. Morita and 
Tatekawa jSHj and Tatekawa et al. jSHj solved the Lagrangian perturbation equations 
for a polytropic fluid up to the second order. Recently, Tatekawa [27^ solved these same 
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equations up to the third order. Hereafter, we call this model the "pressure model." 
Buchert and Dommguez [221 this model the "Euler- Jeans-Newton" (EJN) model. 

In this paper, we study the occurrence of non-Gaussianity in the one-point 
probability density functions (PDFs). For one-point PDFs of the smoothed density 
and velocity fields in a cosmological model, Kofman et al. j2HI analyzed the evolution 
beginning with a Gaussian initial fluctuation. Their analytic results are based on ZA. 
They found that the PDF of the peculiar velocity, both Eulerian and Lagrangian, 
remains Gaussian under linear approximation. For the density fluctuation, they showed 
that the PDF develops a shape similar to a lognormal distribution. 

We consider PDFs as follows. From the point at which the pressure is derived 
from the coUisionless Boltzmann equation (for example, 1201 1^ 1221)5 even if we 
consider only linear perturbation, the pressure model would include nonlinear effects. 
Therefore, even if we analyze the evolution of the peculiar velocity PDF only with linear 
approximation, it will deviate from Gaussian. Please note that the contribution of the 
pressure can be distinguished with PDFs. The effect of the pressure changes the shape 
and the evolution of the PDFs. 

We analyze non-Gaussianity of the PDF of the peculiar velocity and density 
fluctuation. For our analyses, we compute skewness and kurtosis, which are statistical 
quantities for non-Gaussianity of the PDF. First, we show the results for the dust model 
with ZA and N-body simulation. Our results coincide with those of previous efforts. 
Then we indicate that the skewness and the kurtosis are useful for the analysis of non- 
Gaussianity. Next we analyze the PDFs for the pressure model. Here we consider 
both the Eulerian and the Lagrangian linear approximations. In the pressure model, 
because the growing rate of the fluctuation depends on the scale, the occurrence of the 
non-Gaussianity was expected. However, under linear approximation, the effect of the 
pressure does not contribute to occurrence of the non-Gaussianity very much. The PDF 
of the peculiar velocity remains largely Gaussian during evolution. For the PDF of the 
density fluctuation, the Eulerian linear approximation does not indicate the occurrence 
of non-Gaussianity. On the other hand, in the Lagrangian linear approximation non- 
Gaussianity does obviously occur. Furthermore, the evolution of the PDF depends on 
the "equation of state." Therefore if we analyze carefully the PDFs obtained from 
observations, we can expect to find a constraint to the "equation of state" for the 
cosmological fluid. 

This paper is organized as follows. In Sec. |21 we briefly present the evolution 
equation for cosmological fluid. For simplicity, we consider only the Einstein-de Sitter 
Universe model. 

In Sec. we analyze non-Gaussianity of the one-point PDFs. Here we introduce 
two statistical quantities, skewness and kurtosis. First, in Sec. 13.11 we analyze the 
dust model. According to past study, the PDF of the density fluctuation develops a 
shape similar to a lognormal distribution. Its non-Gaussianity could be detected with 
the skewness and kurtosis. In Sec. 13. 2| we analyze the pressure model. In the dust 
model, the PDFs in the Eulerian linear approximation remains Gaussian. Also in the 
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pressure model, this tendency was unchanged. The PDF of the density fluctuation in 
the Lagrangian hnear approximation was sensitive to the "equation of state." 
In Sec. m we discuss our results and state conclusions. 



2. The evolution of fluctuation for the cosmological fluid 

Here we briefly introduce the evolution equation for cosmological fluid. In the comoving 
coordinates, the basic equations for cosmological fluid are described as 

^ + iv..wi + .)}=o, (1) 

dvl, al^l^^ 

— + -{v V^)v + -v = -g V.P , (2 

at a a a ap 

V,xg = 0, (3) 

^x -g = - ATiGpbaS , (4) 

5 (5) 

Pb 

In the Eulerian perturbation theory, the density fluctuation 6 is regarded as a 
perturbative quantity. In linear approximation, from equation (^, we obtain 

^ + -V^-v = 0. (6) 

ot a 

Then we take a divergence of equation Q and substitute equation (jH)) to it. Finally 
we obtain the evolution equation for the density fluctuation in the Eulerian linear 
approximation ^U] : 

§ + 2^|-4.Gp.*-ivJP = 0. (7) 

ot^ a at pbO 

When we assume a polytropic fluid (P = np'^), equation (|7j) becomes 

^^^2^-%-4.Gp,6-^VlS = 0. (8) 

ot-^ a ot 

On the other hand, in the Lagrangian perturbation theory, the displacement from 
homogeneous distribution is considered. 

X = q + s{q,t) , (9) 

where x and q are the comoving Eulerian coordinates and the Lagrangian coordinates, 
respectively, s is the displacement vector that is regarded as a perturbative quantity. 
From equation ©, we can solve the continuous equation (jH) exactly. Then the density 
fluctuation is given in the formally exact form. 

6=l-J-\ J = det(^] . (10) 



\9qjy 

J means the Jacobian of the coordinate transformation from Eulerian x to Lagrangian 
q. Therefore when we derive the solutions for s, we can know the evolution of the 
density fluctuation. 
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The peculiar velocity is given by 

V = as . (11) 

To solve the perturbative equations, we decompose the Lagrangian perturbation to the 
longitudinal and the transverse modes: 

s =VS + s^, (12) 
V-s^ = 0, (13) 

where V means the Lagrangian spacial derivative. Hereafter we consider only the 
longitudinal mode. 

Using the Lagrangian displacement, we obtain the first-order perturbative equation. 

V2 ( + 2^5(1) - AirGp^S^^^ - '^lIpLI^^s^A = . (14) 
\ a a^ J 

Let us compare equations (jH)) and (fT^ . In both of the equations, the same operator 
acts on the perturbative quantity |25j- Therefore in the linear approximation, the form 
of the Lagrangian solutions is the same as that of the Eulerian solutions. 

The first-order solutions for the longitudinal mode depend on spacial scale. 
Therefore the solutions are described with a Lagrangian wavenumber K. For simplicity, 
in this paper, we discuss only perturbative solutions in the Einstein-de Sitter Universe 
model 



S^^\K,t) = C^{K)D+{K,t) + C-{K)D-{K,t) , (15) 
t) - I t-'^'J^^m-e,M\K\t-^-^/^) for 7^ |, 



where JT" means Bessel function. C^{K) is given by the initial condition. 

In this model, the behavior of the solutions strongly depends on the relation 
between the scale of fluctuation and the Jeans scale. Here we define the comoving 
Jeans wavenumber as 

K.^i^f. (IT, 

It depends on the ratio between the comoving Jeans wavenumbers and the wavenumber 
of the fluctuation whether the fluctuation grows. When we ignore the pressure term, 
we obtain the solutions of ZA. 

S(^) =t^/'S+iq)+t-'S_{q). (18) 

S±{q) is given by the initial condition. 
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3. Non-Gaussianity of the one-point distribution functions 

We analyze the non-Gaussianity of the PDF of the peculiar velocity and the density 
fluctuation. To investigate non-Gaussianity, we compute two statistical quantities j3 

EH EH!. 

skewness : 




kurtosis : ( ) ~ 3 , 



cTp = {{p — (p))^) ,P'- physical quantity . 

The skewness and the kurtosis shows display asymmetry and a non-Gaussian degree 
of "peakiness." If the distribution is completely Gaussian, both the skewness and the 
kurtosis become zero. 

For simplicity, we set the Gaussian density field with a scale-free spectrum: 

V{k)(xk. (19) 

To avoid divergence of the density fiuctuation, we introduce a top-hat cutoff at 
k'^ = Ih-^Mpc (at a = 1) in comoving Eulerian coordinates. 

The Gaussian initial condition is generated by COSMICS jSS]- We set up the 
initial condition at a = 10~^. The initial peculiar velocity and the density fiuctuation 
are adjusted by the growing solution of ZA. Both the skewness and the kurtosis of the 
initial condition are less than 10~^. 

For N-body simulation, we execute P^M code. The parameters of simulation were 
given as follows: 

Number of particles : = 128'^ , 

Box size : L = 128/i^^Mpc (at a = 1) , 

Softening Length : e = 0.05/i^^Mpc (at a = 1) . 

For perturbative models, we set the parameters as follows: 

Number of grids : = 128^ , 

Box size : L = 128/i"^Mpc (at a = 1) . 

The dispersion, skewness and kurtosis is computed by numerical method. For 
analyses of non-Gaussianity, we coarsen the fields. For the density fiuctuation, we 
coarsen the density field over Ih-^Mpc (at a = 1) in comoving Eulerian coordinates with 
a top-hat window function. For the peculiar velocity, we analyze on grids at intervals 
of l/i^^Mpc (at a = 1) in Lagrangian space for Lagrangian perturbation models. In 
N-body simulation, we analyze the peculiar velocity of each particle. 

We compute the dispersion, skewness and kurtosis at ten time slices from a = 0.1 
to a = 1.0. The time interval is given as Aa = 0.1. 
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3.1. Non-Gaussianity in the dust model 

First, we analyze the dust model. Here we analyze the PDF of the peculiar velocity 
and the density fluctuation for the Eulerian linear approximation, ZA and N-body 
simulation. For the dust models, Kofman et al. j2H] analyzed one point PDFs in detail. 
They proved that for the PDF of the peculiar velocity, the Lagrangian PDF is equal to 
the Eulerain PDF at all times. For the PDF of the density fluctuation, Padmanabhan 
and Subramanian [HII derived the distribution function in the nonlinear regime using 
ZA. 

Now we analyze the quasi-nonlinear stage. Figure ^ shows the dispersion of the 
density fluctuation. In our model, the dispersion becomes a > 0.1 at a = 0.1. Because 
of shell- crossing, the growth of the density fluctuation becomes slow gradually. 

Figure El shows the PDF of the density fluctuation in N-body simulation. During 
evolution, the PDF approaches log-normal distribution. The result coincides with that 
of the past analyses j2Hl Ell IH2j- Figure El shows the PDF of the peculiar velocity. 
Because we set up almost isotropic initial condition, even if we pay attention to the one 
direction, we think that the generality of the result is not lost. Therefore, we choose 
the x-direction peculiar velocity. During evolution, because the dense region attracts 
surrounding matter, the probability of fast velocity increases. At a = 1.0, both the 
density fluctuation and the peculiar velocity obviously show non-Gaussian distribution. 
Therefore the skewness and the kurtosis must have non-zero value. 

Figure IH shows the skewness and the kurtosis of the density fluctuation. In 
N-body simulation, the PDF of the density fluctuation is well-fitted by log-normal 
distribution [2B|. Therefore N-body simulation realizes strongly nonlinear evolution, 
then the skewness and the kurtosis increases rapidly. In ZA, the skewness and the 
kurtosis also increases during evolution. However, ZA can describe only quasi-nonlinear 
evolution, and after shell-crossing, i.e., the formation of the caustic, ZA cannot describe 
the evolution of nonlinear structure. Therefore the evolution of the skewness and the 
kurtosis stops gently. 

FigureElshows the skewness and the kurtosis of the peculiar velocity. Because we set 
initial condition by almost isotropic distribution, the initial PDF of one direction of the 
peculiar velocity is Gaussian. As a result the skewness and the kurtosis increase rapidly. 
In ZA, the skewness and the kurtosis also increase during evolution. However, ZA can 
describe only quasi- nonlinear evolution, and after shell-crossing, i.e., the formation of the 
caustic, ZA cannot describe the evolution of nonlinear structure. Therefore the evolution 
of the skewness and the kurtosis stops gently. In the Eulerian linear approximation, the 
PDF of the peculiar velocity and the density fluctuation remains Gaussian at all times. 

From our analyses of the PDF in the dust model, we can conclude that both the 
skewness and the kurtosis are useful quantities for analyses of non-Gaussianity. Next 
we adopt these analyses for the pressure model. 
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Figure 1. The dispersion of the density fluctuation in ZA and N-body simulation. 
Because of shell-crossing, the growth of the density fluctuation becomes slow gradually. 
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Figure 2. The PDF of the density fluctuation in N-body simulation. During evolution, 
the PDF approaches lognormal distribution. 
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3.2. Non-Gaussianity in the pressure model 

In this subsection, we analyze the pressure model. Here we choose the polytropic index 
7 = 4/3,5/3. In a previous paper we compared the density field between N-body 
simulation and the Lagrangian approximations. In this comparison, we showed that 
if a small value was chose for the initial Jeans wavelength Kj (equation (fT7|) ). it is 
hard for a nonlinear structure to form the pressure suppress the evolution of the density 
fiuctuation. On the other hand, if we choose a large Kj value, the evolution of the 
density fiuctuation becomes almost same as that in the dust model. Therefore if we 
want to consider the effect of the pressure and the formation of the nonlinear structure. 
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Figure 3. The PDF of the pecuhar velocity in N-body simulation. During 
evolution, because of attraction to high density region, the PDF deviates from Gaussian 
distribution. 
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Figure 4. The skewness and the kurtosis of the density fluctuation. After quasi- 
nonlinear stage (<5niax > l)j the positive fluctuation evolves quickly. Then the 
skewness increases remarkably. Of course N-body simulation realizes strongly nonlinear 
evolution; the skewness and the kurtosis increase rapidly. For ZA, because of shell- 
crossing, the evolution of the skewness and the kurtosis stops gently. 

we should choose an appropriate value for Kj. In this paper, following the results of our 
previous paper [34j, we choose that the initial Jeans wavelength Kj = 64 at a = 10^'^. 
Here we do not specify a dark matter model; we just analyze the behavior of the pressure 
models. In the pressure model, the evolution of the perturbation depends on the scale. 
Therefore even if the initial distribution is Gaussian, we expect that the non-Gaussianity 
would appear during evolution. 

Figure El shows the dispersion of the density fluctuation. In past papers [SHIEHI, we 
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Figure 5. The skewness and the kurtosis of the pecuhar velocity. Here we choose 
one direction of the peculiar velocity. In N-body simulation, the high-density region 
attracts surrounding matter, and the kurtosis increases. On the other hand, when 
we set an isotropic initial condition, the skewness almost does not increase. In ZA, 
because the evolution of the peculiar velocity is always linear, the PDF of the peculiar 
velocity remains Gaussian at all times. 
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Figure 6. The dispersion of the density fluctuation in ZA and the pressure models. 
During evolution, the dispersion in the case of 7 = 5/3 approaches to that in ZA. 



show that the asymptotic behavior of the perturbative solution in the case of 7 = 5/3 is 
similar to that in ZA. Here we show that the statistical quantities in the case of 7 = 5/3 
approach to these in ZA during evolution. As was done with the dust model, we analyze 
the quasi-nonlinear stage in the pressure model. 

Figures [7| and IHl shows the skewness and the kurtosis of the peculiar velocity in the 
pressure model, respectively. Under the linear perturbation, for the peculiar velocity, 
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Figure 7. The skewness of the pecuhar velocity in the pressure model. Here we 
choose one direction of the peculiar velocity. In the dust model, because the evolution 
of the peculiar velocity is always linear, the skewness takes a constant value. In the 
pressure model, although the non-Gaussianity appears during evolution, its extent 
remains small in the quasi-nonlinear stage. 
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Figure 8. The kurtosis of the peculiar velocity in the pressure model. Here we choose 
one direction of the peculiar velocity. As in the case of the skewness, in the dust 
model, because the evolution of the peculiar velocity is always linear, the kurtosis 
takes a constant value. In the pressure model, although the non-Gaussianity appears 
during evolution, its extent remains small in the quasi-nonlinear stage. 



the Lagrangian PDF is equal to the Eulerain PDF at all times [2H1- In the pressure 
model, although the non-Gaussianity appears during evolution, its extent remains small 
in the quasi-nonlinear stage. In the dust model, the PDF of the peculiar velocity never 
changes. Even if the nonlinear structure forms, the motion of the matter does not stop. 
On the other hand, because the pressure especially affects the dense region, the motion 
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Figure 9. The skewness and the kurtosis of the density fluctuation in the Eulerian 
models. Although we consider the effect of the pressure, the PDF of the density 
fluctuation is almost Gaussian during evolution. 



of the matter is slowed down. Therefore, as we show in Figure the skewness of the 
pecuhar velocity in the pressure model grows more than ten times bigger than that in 
ZA. 

Next we show the skewness and the kurtosis of the density fluctuation in the 
pressure model. First we analyze the Eulerian linear approximation. Even if we consider 
the Eulerian linear approximation, the evolution of the fluctuation is not linear, because 
of the pressure. The evolution depends on the scale of the fluctuation. Therefore the 
PDF of the density fluctuation may deviate from Gaussian. Figure El shows the result. 
Although the skewness and the kurtosis oscillate, the value is still small in a nonlinear 
regime. Therefore we conclude that the PDF of the density fluctuation almost retains 
Gaussianity. 

We analyze the PDFs in the Lagrangian linear approximation. As we show 
in equation (jlOj) . the relation between the density fluctuation and the Lagrangian 
displacement is nonlinear. Therefore as we show in the case of the dust model, even if 
we consider the linear approximation, the PDF of the density fluctuation deviates from 
Gaussian. 

Figures Cni and [TTl show the skewness and the kurtosis of the density fluctuation. 
In the dust model, because of shell-crossing, the growth of the skewness stops gently. In 
the pressure model, although the growth of the skewness lasts longer than in the case 
of the dust model, it finally stops because of shell-crossing. The behavior of the growth 
of the non-Gaussianity depends on the "equation of state." 
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Figure 10. The skewness of the density fluctuation in the Lagrangian models. In the 
dust model, because of shell-crossing, the growth of the skewness stops gently. In the 
pressure model, although the growth of the skewness lasts longer than in the case of 
the dust model, it finally stops because of shell- crossing. 
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Figure 11. The kurtosis of the density fluctuation in the Lagrangian models. As in 
the case of the skewness, in the dust model, because of shell-crossing, the growth of the 
kurtosis stops gently. In the pressure model, although the growth of the kurtosis lasts 
longer than in the case of the dust model, it finally stops because of shell-crossing. 



4. Summary 

We analyzed a one-point PDF for the peculiar velocity and the density fluctuation. 
Here we start from the Gaussian distribution and analyze the occurrence of non- 
Gaussianity. For the dust model, in Eulerian linear approximation, because the growth 
of the fluctuation is always linear, the distribution always remains Gaussian. In 
Lagrangian linear approximation, although the PDF of the density fluctuation becomes 
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non-Gaussian, the PDF of the pecuhar velocity remains Gaussian, because the pecuhar 
velocity is proportional to the Lagrangian displacement. Therefore the non-Gaussianity 
is produced by nonlinear effect. 

Next we consider the effect of the pressure. If the pressure exerts an effect, the 
evolution of the fluctuation depends on the effect's scale. Therefore even if we take 
only the Eulerian linear approximation, the pressure may produce non-Gaussianity of 
the density fluctuation. However, from our results, the PDF of the density fluctuation 
in the Eulerian linear approximation is almost Gaussian during evolution. In the same 
way, the PDF of the peculiar velocity is almost Gaussian during evolution. Therefore we 
can conclude that most of the non-Gaussianity of the PDF is produced by the nonlinear 
effect of the perturbation. 

We compared the evolution of the non-Gaussianity in the Lagrangian linear 
approximation. As we showed in Figures HOI and [HJ the evolution depends on the 
"equation of state." In analyses of the PDFs of both the peculiar velocity and the density 
fluctuation, we can distinguish the polytropic index. According to past analyses |2SII2S1, 
the behavior of the linear perturbative solutions differ widely between the cases of 
7 = 4/3 and 7 = 5/3. In the case of 7 = 4/3, the behavior of the solutions strongly 
depends on the scale of the fluctuation. If we consider large-scale fluctuation, because 
self-gravity dominates, the fluctuation grows. On the other hand, if we consider small- 
scale fluctuation, because the pressure dominates, the fluctuation oscillates and decays. 
In the case of 7 = 5/3, the pressure affects only the early stage. At the late time, 
the behavior of the perturbative solution approaches that of ZA solutions. Therefore 
the statistical quantities we computed, i.e., the variance, the skewness and the kurtosis 
of the peculiar velocity and the density fluctuation approaches to those in the case of 
ZA (Figures H El H CUl and HH) . 

Although our result can apply until quasi-nonlinear stage is reached, we can expect 
that an "equation of state" can be distinguished from the growth of the non-Gaussianity. 
If we can observe the PDF of the density fluctuation in high-z region, i.e., the quasi- 
nonlinear region, we can flnd the constraint of the character of the dark matter [22] • 

For comparing between the theoretical models and the observation, spacial two- 
point correlation function jHHl IHZl IHHl IHEl HOI HH Il2j can also be important. We will 
compute the correlation function and discuss the difference between the theoretical 
models. 
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